Quantum two-level systems in Josephson junctions as naturally formed qubits 
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The two-level systems (TLSs) naturally occurring in Josephson junctions constitute a major ob- 
stacle for the operation of superconducting phase qubits. Since these TLSs can possess remarkably 
long decoherence times, we show that such TLSs can themselves be used as qubits, allowing for a 
well controlled initialization, universal sets of quantum gates, and readout. Thus, a single current- 
biased Josephson junction (CBJJ) can be considered as a multiqubit register. It can be coupled to 
other CBJJs to allow the application of quantum gates to an arbitrary pair of qubits in the system. 
Our results indicate an alternative way to realize superconducting quantum information processing. 



Several advances in the field of quantum informa- 
tion processing using superconducting circuits have been 
made in recent years 0. A major obstacle to further ad- 
vances, however, is the problem of decoherence. In par- 
ticular, recent experiments on current-biased Josephson 
junctions (CBJJs) revealed resonances that suggest the 
presence of quantum two-level systems (TLSs) that are 
strongly coupled to the CBJJ when it is biased near one of 
those resonances 0,0- So far the TLSs have been treated 
as a nuisance that prevent the operation of the CBJJ as a 
qubit near any resonance. In this paper we show that the 
TLSs themselves can be used as qubits. The CBJJ then 
acts as a bus, enabling state initialization, one- and two- 
qubit operations, and readout. Moreover, the results of 
Refs. [2 Q show that the decoherence times of the TLSs 
are longer than those of the CBJJ. That property can 
be used in a scalable design such that the decoherence 
time of the entire system scales as the CBJJ decoherence 
time T^ BJJ , as opposed to the usual T ( f B ' 1 J /N , where N 
is the number of CBJJs in the circuit. Our results there- 
fore demonstrate an alternative way to achieve a scalable 
qubit network in a superconducting system, in addition 
to illustrating a method to perform multi-qubit experi- 
ments with available experimental capabilities. 

Model. — The phase qubit, which is comprised of a sin- 
gle CBJJ, is one of the simplest experimental implemen- 
tations of a qubit in superconducting systems 0,3- The 
working states, |0) and |1), are the (metastable) ground 
and first excited states in a local minimum of the wash- 
board potential produced by the CBJJ. The (undriven) 
system is described by the Hamiltonian 
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Here Q is the operator of the charge on the junction, 
C is the junction's capacitance, (p is the operator of 
the Josephson phase difference, and I c ,h(< Ic) are 
the critical and bias current, respectively. The non- 
linearity of the potential allows to consider only the 
two lowest energy states. The Pauli matrix o~ z oper- 
ates in the subspace {|0), |1)}. The transition frequency 



wio ~ u p = (27t/ c /$ C) 1/2 (1 - j 2 ) 1/4 , the plasma fre- 
quency in the biased junction (the corrections due to 
nonlinearity are ~ 10% and j = Ib/I c (see e.g. Q). 
The terms in the Hamiltonian proportional to a x and 
a y , which enable a complete set of onc-qubit gates, ap- 
pear in the rotated frame of reference when applying mi- 
crowave pulses of bias current at the resonance frequency, 
h -* foe + Iftw^t)cos{u}\Qt) + Iftwstysinfaiot), as ex- 
plained in Ref. |5[. 

The simplicity of the qubit design and manipula- 
tion contributed to its successful experimental realiza- 
tion [H Q and a spectroscopic demonstration of the 
formation of entangled two-qubit states Q . Nevertheless 
the experiment in Ref. Q also demonstrated that TLSs 
present in the tunneling barrier tend to destroy the co- 
herent operation of the qubit (such as Rabi oscillations). 

Such TLSs are ubiquitous in solid state systems wher- 
ever disorder is present and can be thought of as groups 
of atoms capable of tunneling through a potential bar- 
rier between two degenerate configurations. They are 
currently believed to be the main source of 1/ /-noise in 
solids. 

The observed coherent oscillations between a TLS and 
a phase qubit proved that a TLS can be considered as 
a coherent quantum object described by the pseudospin 
Hamiltonian 

Htls = ~ ycr* - |cr z , (2) 

with a decoherence time longer than that of the qubit. 
The Pauli matrices d x and a z operate on the TLS states. 
Note that, since the nature of the TLSs is currently un- 
known, one cannot derive the values of the TLS param- 
eters from first principles. As will become clear below, 
however, neither a derivation of those parameters nor an 
understanding of their physical origin is necessary in or- 
der to make use of the TLSs. 

Although the available experimental data gives rather 
limited information about the TLS-TLS interaction, it is 
highly unlikely that the TLSs interact directly with each 
other or with the external fields, because of their sup- 
posedly microscopic dipole moments and relatively large 
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spatial separation. 

The TLS-CBJJ coupling is believed to be due to one 
of the following mechanisms: (A) through the critical 
current dependence on the TLS position |2| or (B) the 
direct dipole coupling to the junction charge Q (which is 
currently considered more likely) 0,0: 
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Both produce the coupling term h x a x a z . In addition, 
in the former case there appears an h z a z a z -teiui, which 
reflects the change in the inter-level spacing u>i due 
to the plasma frequency dependence on the bias. The 
ratio 7 ee h z /h x = {E c l2Ejfl\\ - where 
E c = 2e 2 /C and Ej = I c <£ /2ir. If (E c /2Ej) = 10" 7 
(achievable using the external capacitor technique Q) 
and j = 0.90, then 7 = 0.05 Q. 

Hereafter, we consider the case where 7 < 1. Oth- 
erwise, the effective TLS-CBJJ decoupling is impossible 
even when out of resonance (on the positive side, this 
effect would provide the means to definitively establish 
the mechanism of TLS-CBJJ coupling). Then, perform- 
ing a basis transformation on the TLS states, relabeling 
[A 2 + e 2 ] 1 / 2 — > A, and using the rotating- wave approx- 
imation, we finally obtain the effective Hamiltonian for 
the TLS-CBJJ system: 
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with the effective coupling coefficients Xj. The coupling 
term acts only in resonance, when \htuio — Aj \ < Xj. 

The Hamiltonian Q was derived under the experimen- 
tally relevant assumption that the coupling A <C frujio ~ 
A. We neglected a term of the form [ia x a z in Eq. be- 
cause its influence on the system's dynamics is negligible 
when \i <C hwio and n 2 /Tiujiq <C A. 

Qubit operations. — The inter-level spacing hu>io is 
tuned by the bias current, which together with the res- 
onant behaviour of the TLS-CBJJ coupling allows the 
independent manipulation of each TLS (due to the nat- 
ural dispersion of their characteristic energies, Aj). The 
condition A <C A, hujio allows us to consider changes in 
the bias current as adiabatic from the point of view of in- 
ternal CBJJ and TLS evolution, but instantaneous with 
respect to the CBJJ-TLS dynamics. 

First, consider the singlc-qubit operations, assuming 
for the time being that the decoherence times of the 
TLSs, Tp\ exceed the decoherence time of the CBJJ, 
Td BJJ ' which in turn is much larger than the characteris- 
tic interaction time H/X. On resonance with the jth TLS, 
the Hamiltonian ||3J contains a block that acts as — Xja x 
in the subspacc of degenerate states {|1) ® \g), |0) <g> \e)}. 
Its operation leads to quantum beats between the states 



of the TLS and the CBJJ (see Fi; 
riod Tj ee fi/Xj, as was observed 
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qubit operations on a TLS and its initialization to an 
arbitrary state can be achieved as follows: After initial- 
izing the CBJJ in the state |0) we bring it in resonance 
with the TLS for the duration Tj/2; as a result the states 
of TLS and qubit are swapped, 
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Then the CBJJ is taken out of resonance with the TLS, a 
rotation of its state is performed, and the resulting state 
is again transferred to the TLS (always compensating for 
the parasite phase shifts f (1 + A/A)). The readout of the 
qubit state can be performed after the swap ©, e.g. by 
using the technique of Q, where the transition to the 
resistive state, with its potentially-problematic coupling 
to the TLSs, is avoided. The decoherence time can 
be determined using a similar sequence of operations, by 
initializing the TLS in a superposition state, decoupling 
it from the CBJJ, and measuring the probability to find 
the TLS in a given state after a given time. Note also 
that if the Rabi frequency of CBJJ greatly exceeds the 
TLS-CBJJ coupling, the manipulations with the state of 
the CBJJ can be performed while in resonance with a 
TLS, reducing the overall operation time [Icj . 

Note that as soon as the CBJJ is biased away from 
the resonance with the TLS, they effectively decouple, 
and any perturbation of the quantum state of the CBJJ 
does not affect the survival of quantum coherence in the 
TLS. Therefore the quality factor (the number of gates, 
of average duration T gat o, that can be performed in the 
system before it loses quantum coherence) is no less than 
^c? BJJ / T gatc, and can exceed this value depending on the 
specific decoherence decay law. 

The application of two-qubit gates between two TLSs 
inside the same CBJJ can be achieved similarly. For 
example, the universal ISWAP gate can be performed 
by the following sequence of operations: tune the CBJJ 
in resonance with TLSi for Ti/2, then tune it in reso- 
nance with TLS2 for T2/2, and finally again in resonance 
with TLSi for ri/2. (In the process each TLS qubit was 
also rotated by tt/2 around z-axis, which can be trivially 
repaired by performing the single-qubit rotation as de- 
scribed in the previous paragraph.) At the end of the 
above sequence, the CBJJ is decoupled and disentangled 
from both TLSs. 

Numerical estimates. — The idealized picture of two- 
qubit operations neglected (1) the finite detuning be- 
tween TLS qubits, (2) the influence of other TLSs present 
in the CBJJ, (3) decoherence and (4) the finite time of 
the CBJJ bias current adjustment. We numerically in- 
vestigated the impact of these effects on the fidelity F 
defined as the minimum 
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where is an initial state. C/ideal is the (ideal) desired 
operation, and p/ is the numerically obtained final den- 
sity matrix, (note that the quantity inside the parenthese 
in Eq. JS} depends on the initial state). 
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FIG. 1: Fidelity of the ISWAP gate on two TLS qubits in 
a CBJJ as a function of the relative difference between the 
TLS energy splittings (Ai and A2, respectively). The solid 
line corresponds to no decoherence, the dashed and dotted 
lines correspond to CBJJ decoherence rates Fi = A/207T and 
Ti = \/2ir, respectively, where in both cases F2 = 2Fi. The 
much slower decoherence of the TLSs have been neglected. 
A = 0.005Ai. (Inset) Schematic depiction of the quantum 
beats between the CBJJ and the TLSi in resonance; TLS2 is 
effectively decoupled from the CBJJ. 

In the numerical simulations, we took the ISWAP op- 
eration as a representative quantum gate and used 900 
different initial states. The results are as follows (see 
Fig. HJ). For two TLS qubits, with Ai = A 2 = 0.005Ai 
(values close to the experimental data the fidelity 

first reaches 90% when S = |Ai - A 2 | = 3.5A, and 99% 
when 6 w 10A. When 5 = 8A, the fidelity is 98.8%. 
Adding an idle TLS with A = 0.002Ai in resonance with 
one of the qubit TLSs reduces the fidelity to 80% (94% 
for A = O.OOlAi). Finally, with two qubit TLSs with 
reduced energies 1 and 1.04, four idle TLSs with reduced 
energies 0.99, 1.01, 1.03 and 1.05, and the reduced cou- 
pling strengths A qu bit = 0.005 and Aidie = 0.002, respec- 
tively, we find that the fidelity of the ISWAP gate « 95%. 
Therefore, we conclude that the finite detuning and pres- 
ence of idle TLSs per se is not dangerous. 

We now take into account the decoherence in the 
CBJJ, neglecting the much weaker one in the TLSs. The 
ISWAP fidelity with 5 = 8Xj and without idle TLSs is 
81% when Ti = (l/2)r 2 = A/20tt, but it drops to 13% 
when Ti = (l/2)r 2 = A/2tt. 

Finally, we now consider the effects of a finite CBJJ 
bias switching time between the resonant frequencies of 



the TLSs. Taking a simple linear ramp with t — 2 ns 
and no pulse optimization, and neglecting decoherence, 
we find that the fidelity drops to 80%. Finite decoher- 
ence Ti = (l/2)r 2 = A/20-7T further suppresses it to 63%. 
The above estimates show that the operation of the pro- 
posed two-qubit gate can be realized with the current 
experimental techniques used, e.g., in Ref. 

Scalability. — To ensure the scalability of the system, 
we must be able to perform two-qubit gates on TLSs lo- 
cated in different CBJJs. It can be done by swapping the 
states of the capacitively coupled adjacent CBJJs like in 
[i~l| . but a better solution is based on the method sug- 
gested in • Here the qubit-carrying CBJJs are coupled 
capacitively to a common linear LC circuit with reso- 
nance frequency luq much higher than the characteristic 
frequencies of separate CBJJs (Fig. EJ. 

Using the standard formalism, we can write the La- 
grangian of the system, 
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where <j)j (t) = Vj (t) , the voltage between node j and the 
ground. The corresponding Hamiltonian becomes 



1 N 1 N LU 2 



j,k=0 



2 ^ C~ 

3=0 ^33 



where (j> = (h/2e)(p (see Eq. JQ), [<f>j,Pj] = ifi, C _1 is the 
inverse capacitance matrix, ujq is the frequency of the LC 
bus, ujj is the frequency of the jth CBJJ in the harmonic 
approximation, and the ellipsis stands for the nonlinear 
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FIG. 2: Scalability of the structure: Two-qubit operations 
between the TLSs on different CBJJs (j = 1, . . . ,N) are en- 
abled by the common LC circuit, which is capacitively coupled 
to the CBJJs. The Josephson energy, capacitance, and bias 
current of the jth CBJJ are Ej, Cj, Ij resp. 
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corrections. After introducing the Bose operators a and 
a) via tfij = A.j(a,j + a}-)/2,f>j = h(a,j — at)/(iAj), Aj = 

1/2 

({2hC~j 1 )/ujj) , the Hamiltonian becomes 
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with the effective coupling 
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By assumption, too 3> Wj, j ' = 1, . . . AT. Therefore the 
Hamiltonian can be projected on the ground state of the 
LC bus • The nonlinearity of the CB J J allows us to 
further restrict the Hamiltonian to the subspace spanned 
by the states |0), |1) of each CBJJ, producing 



H, 



eff — - 



1 N 



N 

E 

k>j=i 



(8) 



In the interaction representation with respect to Hq = 
\ Ylj= i a i it is easy to see that a y (t) = a y (0) cos Uj t+ 

(0) sinwji 0. Therefore, the pairwise couplings in (JHJ 
will only be effective for the CBJJ tuned in resonance 
with each other, in which case (ujj = u>k) the effective 
interaction term is H^g = gjk i"',^', + ^l^y) /2- On the 
subspace spanned by the states \0)j ® |l)k and ® 
|0}fe the operator (o-jcr^ + &yGy) /2 acts as (while it 
is exactly zero outside). Therefore this coupling allows 
the same universal two-CBJJ manipulations as in [ill] . 
Universal two-qubit gates on TLSs situated in different 
CBJJs can then be performed by transferring the states 
of the TLS1.2 to the corresponding CBJJ1.2, performing 
the two-qubit operations on the states of CBJJ12, and 
retransferring the resulting states back to the TLSs. 

The number of TLS per CBJJ is of order 10 and de- 
pends on the fabrication. The decoherence times of TLSs 
are determined by their local environment and are insen- 
sitive to the number of CBJJs linked to the same LC cir- 
cuit. Similarly, the influence of the LC circuit on the de- 
coherence of the CBJJ is negligible as long as luq 3> u>io ~ 
uj p , the same requirement we need to obtain the cou- 
pling H?.ff- Therefore the scalability of the system is lim- 
ited by the condition [L(C + NC[ . off )] _1 > 7 c /[*oCi], 
or, to the same accuracy, N <C {&o/LI c ] [Ci/C[] . There- 
fore N can be of order few dozen without violating the 
applicability conditions for the above scalable design. 

The usefulness of TLSs for our purpose could be ques- 
tioned because their parameters undergo spontaneous 
changes. Nevertheless, since such changes typically hap- 
pen on the scale of days |l5j . a pre-run calibration is 
sufficient for any realistic task. Another concern that 
could be raised is the fact that any operation on the 
TLSs is done through the CBJJ, so that the latter's de- 
coherence must be a limiting factor on the number of 



operations that can be performed on the TLSs, no mat- 
ter how long their decoherence times are. However, an 
important point to note here is that in the above design 
with N CBJJs, only one or two of them are used during 
any gate operation. Therefore, the decoherence time of 
the entire system is of the order of T^ BJJ rather than 
31 /N, which one would normally obtain when using 
N CBJJs as phase qubits. 

It is worth noting that there is some control, albeit very 
limited, over the properties of the TLSs (e.g. the number 
of TLSs in a CBJJ was drastically reduced by using a 
different barrier material 0). They will not be the only 
naturally formed objects to allow quantum manipulation. 
A case in point is the observation of Rabi oscillations in 
self-assembled quantum dots 0] . 

Summary. — In conclusion, we have demonstrated that 
quantum two-level systems (TLSs) naturally occurring 
in CBJJs can be used as qubits. The one- and two- 
qubit gates, initialization and readout can be readily per- 
formed, and the system can be scaled beyond a single 
CBJJ. Being microscopis objects, TLSs have a higher 
probability of possessing long decoherence times, which 
are only weakly affected by the TLS-CBJJ interactions, 
due to that interaction being switched off for most of 
the time. The tunability of the CBJJs compensates for 
our currently limited control over the parameters of the 
TLSs. In any case, both our numerical simulations and 
especially the observation of quantum beats between TLS 
and CBJJ Q show that the experimental realization of 
our scheme is within the reach of current experiments 0] • 

We are grateful to M. Grajcar and Y. X. Liu for valu- 
able comments. This work was supported in part by 
the ARO, LPS, NSA and ARDA under AFOSR con- 
tract number F49620-02-1-0334; and also supported by 
the NSF grant No. EIA-0130383. A.Z. acknowledges par- 
tial support by the NSERC Discovery Grants Program. 
S.A. was supported by a fellowship from the JSPS. 



[1] 
[2] 

[3] 
[4] 
[5] 
[6] 
[7] 
[8] 

[9] 



[10] 
[11] 



J. Q. You and F. Nori, Phys. Today 58 (11), 42 (2005). 
R.W. Simmonds et al, Phys. Rev. Lett. 93, 077003 
(2004). 

J.M. Martinis et al, Phys. Rev. Lett. 95, 210503 (2005). 
J.M. Martinis et al., Phys. Rev. B 67, 094510 (2003). 
J.M. Martinis et ai, Phys. Rev. Lett. 89, 117901 (2002). 
A.J. Berkley et al., Science 300, 1548 (2003). 
K.B. Cooper et ai, Phys. Rev. Lett. 93, 180401 (2004). 
I. Martin, L. Bulaevskii, and A. Shnirman, Phys. Rev. 
Lett. 95, 127002 (2005). 

Note that the temperature does not have to be smaller 
than E.j (with the proper unit conversion). It only has 
to be smaller than the typical inter-level spacing (~ 100 
mK). 

S. Ashhab, J.R. Johansson, and F. Nori, New J. Phys. 8, 
103 (2006). 

A. Blais, A. Maassen van den Brink, and A.M. Zagoskin, 



5 



Phys. Rev. Lett. 90, 127901 (2003). 1999), p.301. 

[12] Yu. Makhlin, G. Schon, and A. Shnirman, Nature 398, [14] A. Miiller et al, Appl. Phys. Lett. 84, 981 (2004). 

305 (1999). [15] R. McDermott, private communication. 
[13] See e.g. M. Orszag, "Quantum Optics" (Springer, Berlin, 



